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ABSTRACT 

We present an integral field study of the internal structure, kinematics and stellar population 
of the almost edge-on, intermediate luminosity {L^) elliptical galaxy NGC 4697. We build 
extended 2-dimensional (2D) maps of the stellar kinematics and line-strengths of the galaxy 
up to ^ 0.7 effective radii (i?eff) using a mosaic of 8 VIMOS (Visible Multi-Objects Spectro¬ 
graph on the VLT) integral-field unit pointings. We find clear evidence for a rotation-supported 
structure along the major axis from the 2D kinematical maps, confirming the previous clas¬ 
sification of this system as a ‘fast-rotator’. We study the correlations between the third and 
fourth Gauss-Hermite moments of the line-of-sight velocity distribution (LOSVD) and 
with the rotation parameter (V/a), and compare our findings to hydrodynamical simulations. 
We find remarkable similarities to predictions from gas-rich mergers. Based on photometry, 
we perform a bulge/disk decomposition and study the stellar population properties of the two 
components. The bulge and the disk show different stellar populations, with the stars in the 
bulge being older ( agebuige = 13.51 Gyr, agedisk = 10.5l2‘.o^yr) and more metal-poor ( 
[M/H]buige = — O.lTlo.i^^ [M/H]disk = — 0.03lo!i^)- The evidence of a later-formed, more 
metal-rich disk embedded in an older, more metal-poor bulge, together with the LOSVD struc¬ 
ture, supports a mass assembly scenario dominated by gas-rich minor mergers and possibly 
with a late gas-rich major merger that left a previously rapidly rotating system unchanged. 
The bulge and the disk do not show signs of different stellar Initial Mass Function slopes, and 
both match well with a Milky Way-like IMF. 

Key words: galaxies : formation - galaxies : evolution - galaxies : elliptical and lenticular 


1 INTRODUCTION 


Integral Field Spectroscopy (IFS) makes it possible to derive spa¬ 
tially resolved and continuously sampled line-of-sight velocity dis¬ 
tribution (LOSVD) 2-dimensional (2D) maps, providing a unique 
view of the internal structure of galaxies, as well as detailed mor¬ 
phological and dynamical information, which are essential to fully 
understand the physical mechanisms that regulate galaxy formation 
and evolution. 

An increasing number of studies have recently derived the in¬ 
tegrated and spatially resolved properties of a statistical sample of 
nearby galaxies using IFS (e.g. Bacon et al.||2001[ [Thatte et al. 


2006 JRawle et al.|2008[|Czoske ct al.|2008| Scott et al.|2012[[Fog¬ 
arty et al.|2014[|Gonzalez Delgado et al.|2014||Sanchez-Blazquez 
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|et al.|2014 [Walcher et al.|2014ll. Among the most successful, the 
ATLAS Cappellari et al.]207l] l data set provides a complete in¬ 
ventory of the baryon budget and a detailed 2D analysis of stel¬ 
lar and gaseous kinematics, together with resolved stellar popula¬ 
tions of a complete and statistically significant sample of early-type 
galaxies (ETGs) in the local Universe. 

One of the most interesting results of the ATLAS^^ team 
is that it is possible to connect the 2D features of the line-of- 
sight-velocity distribution (LOSVD) to galaxy formation scenar¬ 
ios by comparing these to high-resolution simulations (e.g. |Naab| 
|et al.|2014j hereafter N+14). For the majority of early-type galax¬ 
ies (~ 86%), the formation processes result in disk-like objects 
that maintain the axisymmetric shape of the progenitors, while only 
a small fraction of systems (the most massive ones, 12%) rotate 
slowly with no indication of disks. Possible processes that can form 
embedded disk components include minor mergers, major merger 
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with disk regrowth, gas accretion events and secular evolution ( |Kra-| 
[jnovic et al.|2011| l. 

It is possible, in principle, to discriminate between the dif¬ 
ferent formation scenarios by accurately analyzing the higher mo¬ 
ments of the LOSVD, the Gauss-Hermite coefficients and 
(representing skewness and kurtosis, respectively) and comparing 
these to high-resolution simulations. Gas has a strong impact on 
the non-Gaussian shapes of the line-of-sight velocity distributions 
(LOSVDs) of the merger remnants. In particular, when the values 
of hs and Vrot have opposite signs, the leading (pro-grade) wing 
is steeper than the trailing (retrograde) wing. This indicates the 
presence of an embedded disk (see e.g. [Bender, Saglia & Gerhard] 
1994| ). When an anti-correlation between the hs and theV/a values 


does not exist but the galaxy shows a fast rotation profile, the rota¬ 
tion could be have caused by a recent gas poor major merger that 
typically does not permit re-growth of a significant disk compo¬ 
nent ( Barnes & Hernquist|1996[ [Naab, Khochfar & Burkert|2006[ 


[Hoffman et al.||2009| hereaft^H-i-09). The sign of the coeffi¬ 

cient distinguishes a peaked shape (h^ > 0) from a flat-topped one 
(/i 4 < 0), where the peak is broad and the wings are narrow. 

The age and metal content of galaxies can provide additional 
important information about the formation mechanism and its time- 
scale. To obtain an unbiased and detailed understanding of the star 
formation history of ETGs from their integrated light, it is neces¬ 
sary to study carefully their unresolved stellar populations, and to 
break the degeneracies existing between age, metallicity and el¬ 
ement abundances (e.g. |Worthey||1994[ |. The best way to achieve 
this is to look at as many stellar features as possible ( [Spiniello] 
[etal.|201^ hereafter S-i-14). Clearly, in this context, the very narrow 
wavelength coverage of SAURON is not ideal, since stellar popu¬ 
lation studies would rely on only a few well known blue-optical 
absorption lines (H/3, Mg6, Fe5015 and Fe5270). 

The VIMOS IFU (Visible Multi-Objects Spectrograph on the 
VLT) covers a wider wavelength range (5250-7400 A using the 
blue grism and 6450-8600 A using the red grism, in its high- 
resolution mode), with comparable spectral and spatial resolution, 
and is clearly a better choice. This point is particularly relevant 
after the recent finding that the low-mass end of the Initial Mass 
Function (IMF) is not universal and might become steeper with in¬ 


creasing stellar mass and stellar velocity dispersion ([van Dokkum 

& Conroy 12010| [Treu et al.|2010| [Spiniello et al.|2012| C 

appellari 

et al.|2012l[2013(|Dutton et al.[201 3a [Spiniello et al.|2014 

La Bar- 

bera et al.[2013| Tortora, Romanowsky & Napolitano [20 13 

Tortora 

et al.|2014b|a| but see also [Smith, Lucey & Conroy|2015 

or with 


changes in other galaxy structural parameters (e.g. compactness, 
stellar density, size, see Spiniello et al. 2015b, submitted). Indeed, 
to constrain the low-mass end of the IMF, it is crucial to look at the 
red part of the optical spectra, where most of the gravity-sensitive 
absorption features are. Some of the most promising indices, like 
the optical TiO molecular lines (bTiO, aTiO, TiOl; S-i-14), can 
be captured with VIMOS but are outside the SAURON coverage. 
Moreover, the SWELLS survey ( [Treu et al.|2011[ ) has shown that 
there could be one universal normalization of the IMF for disk¬ 
like stellar populations (presumably Chabrier-like) and a heavier 
(Salpeter-like) one for the older and possibly more metal rich stel¬ 
lar populations found in massive bulges and spheroids ([Dutton et al.[ 
[2013b[ l. 

Furthermore, 2D integral-field spectroscopy is crucial for con¬ 
straining the low-mass end of the IMF slope, given the very recent 
observational evidence that seems to indicate that the IMF nor¬ 
malization varies within individual galaxies ( [Martm-Navarro QtH] 
[2015[ Sarzi, Spiniello & Krajnovic, 2015, in prep.) following the 


velocity dispersion profile (but see also [Smith, Lucey & Conroy] 
[2015[ [spiniello et al.|2015[ l. In this context, a strong observational 
limitation is the small field-of-view (FOV) of IFS instruments, typ¬ 
ically of the order of tens of arcsec, preventing good coverage of 
the target galaxies on the sky and a proper spatially-resolved study 
of the IMF slope. Only very recently has a new class of integral- 
field spectrographs with larger field of view been built to achieve 
the challenging task of studying the internal structure of galaxies. 
Unfortunately most of them still cover a small range of wavelength 
and therefore are not suitable for constraining the low -mass end 


of the IMF normalization (e.g. VIRUS-P, 3600-6800 A, Hill et al. 


2008 or CALIFA, 3750-7000 A, Sanchez et al. 2012 in its red con¬ 


figuration). The VLT Multi Unit Spectroscopic Explorer (MUSE) 
represents the perfect solution, combining a large field of view 
(1x1 arcmin) and a wide wavelength range (4650-9300 A). Fi¬ 
nally, in the near future the MaNGA (Mapping Nearby Galaxies 
at APO, [Law & MaNGA Team|2014j l project will further improve 
the situation, covering the wavelength range 3600-10000 A up to 
~ 1.5i^eff with a spatial resolution of ~ 2.5 arcsec (and a reso¬ 
lution of R~ 2000). We refer to Table 1 of [Brodie et al.[ ( [20T4] ) 
for a detailed comparison between current and future integral-field 
surveys. 

Thus, at the present day, VIMOS remains among the most 
competitive instruments for galaxy structure and stellar population 
studies, provided that the small FOV is compensated for by mo¬ 
saic observations, i.e. contiguous IFU pointings astrometrized and 
stacked together to create a larger final FOV. 

This is indeed the approach we decided to follow in order to 
attempt to systematically study the stellar population and the kine¬ 
matics of the whole surface of nearby ETGs up to almost the effec¬ 
tive radius over a large wavelength range and with medium-high 
spectral resolution. Given the challenges in the data reduction, cal¬ 
ibration and visualization, with this paper we present the analy¬ 
sis of one of the galaxies observed in a pilot program focusing 
on the kinematics and the stellar population analysis of two bat¬ 
tened galaxies (NGC 821 and NGC 4697|^ In particular we focus 
on NGC 4697, for which we have 8 VIMOS pointings covering 
~ 80" X 60" effective area, which allows us to investigate regions 
up to ~ 45" from the galaxy center, i.e. ~0.7 i^eff (7?eff= 66” 
from [de Lorenzi et al.|2008] l . 


1.1 Target selection and outline of the paper 

NGC 4697 is an E6, almost edge-on galaxy located along the Virgo 
southern extension, with evidence of substructures (from HST ACS 
images, [Lauer et al.|1995| l and the presence of a strong streaming 
motion component (known since the pioneering work of [Bertol^ 
[& Capaccioli||1975] l. The galaxy is representative of the class of 
disky, fast-rotating elliptical L* galaxies ([Bender, Saglia & Ger¬ 


hard] 19941), which are candidates for having less dark matter (jCc 


one & Napolitano[[2004[ [Napolitano et al. 12005] [Cappellari et al. 

2006| ) than the boxy, slow-rotators of comparable mass. As such, 
it is among the sample of elliptical galaxies whose planetary neb¬ 
ula (PN) kinematics show surprisingly little sign of a massive dark 
matter halo ([Mendez et al.[2001[[Romanowsky et al.|2003[[Napoli-[ 
[tano et al.|2005| [Douglas et al.|2007[ [Mendez et al.|2009[ l. Finally, 

the spheroidal part of NGC 4697 may well be a large bulge ( [De-[ 
[jonghe et al.|1996| l and therefore this system is an interesting can¬ 
didate for testing the SWELLS ( [Treu et al.|2011[ | claim about the 

^ Program IDs P079.B-0622A and P085.B-0949 (P.I. N.R. Napolitano). 
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different normalization of the low-mass end of the IMF slope (uni¬ 
versal, Milky Way-like IMF for the disks and more bottom-heavy 
IMF, possibly steepening with velocity dispersion for bulges, [DuF] 
[ton et aL|2013b| l. 

Based on the combined information of the 2D kinematical 
characterization of the system, and stellar population analysis, we 
are able to fully characterize the stellar population of the disk and 
bulge components of NGC 4697. This is important per se in order 
to investigate the variation of the stellar populations and IMF slope 
in different structural components of the system and possibly to in¬ 
fer a formation scenario for it, but it will also allow us to correctly 
normalize the stellar mass of the galaxy in forthcoming dynamical 
analysis. 

The paper is organized as follows: in Section we outline 
the observations and data reduction, in Sectionwe introduce our 
kinematical technique and we show 2D kinematical maps as well as 
1-dimensional (ID) kinematics profiles along the major axis of the 
galaxy. We compare the latter with kinematics profiles in literature. 
In particular, in this Section, we highlight correlations between the 
/i 3 and /i 4 parameters and V/a and try to connect these to predic¬ 
tions from galaxy formation simulations. In Section 4, we focus on 
the stellar population analysis performed via line-index strengths 
measurements on the ID spectra of the disk and of the bulge sep¬ 
arately. We also attempt to constrain the IMF normalization of the 
two components. Finally, discussion and general findings are pre¬ 
sented in Section 5 together with conclusions. 


2 OBSERVATION AND DATA REDUCTION 

The integral-field spectroscopic (IFS) observations have been car¬ 
ried out with VIMOS, the multi-mode wide-field optical instrument 
mounted on the ESO VLT, as part of a proposal carried on in two 
different runs : May 2007 and March 2010, respectively. 

The Integral-field Unit (IFU) head of VIMOS consists of a 
square array of 80 x 80 micro-lenses with a sampling of 0.67'' 
per spatial element (spaxel) at the Nasmyth focus of the telescope. 
In its high-resolution mode only the central 1600 lenses are used, 
resulting in a field-of-view (FOV) of 27” x 27”. The seeing was 
generally below 0.8" during the observations, so that each fiber 
of the IFU gives a spectrum that is essentially independent of its 
neighbors. All data sets have been taken with the high-resolution 
blue grism, covering a wavelength range between 4100 — 6300 A 
with a dispersion of 0.51 A/pixel. The spectral resolution of the 
HR-Blue grism is A/A A = 2550, where A A is the full width at 
half-maximum (FWHM) of the instrument profile. 

In order to cover the effective radius and to have a high S/N 
among the whole FOV, we use a mosaic approach. During the first 
observation run, we have obtained 6 different pointings, each made 
by six exposures of ~ 800 seconds, while for the second run we 
have obtained 2 pointings along the major axis of the system, with 4 
exposures of 825 seconds each. The Position Angle of all pointings 
has been fixed to PA = —20. For accurate measurements of the 
stellar kinematics, one must properly evaluate the sky background 
contamination. Since VIMOS does not have sky-dedicated fibers, 
and the target galaxy covers the entire field-of-view, it has been 
necessary to make separate telescope pointings (of ~ 300 seconds) 
to measure the sky background for each observing night. 


2.1 Data Reduction 

Data reduction has been carried out using the VIMOS feline 
(version 2.1.1) by the ESO “Data Elow System Grou^’, via 
the Gasgano data file organizer developed by ESO. A schematic 
representation of the data reduction processes performed within 
Gasgano is provided by |Zanichelli et al.| ( |2005| ). The reduction in¬ 
volves standard steps for fiber-based integral-field spectroscopy: (i) 
Subtraction of a master bias frame created by averaging all the bias 
frames (ii) Eiber identification (correct association of a fiber posi¬ 
tion on the lEU mask to a corresponding position on the CCD) and 
correction of the differences in the fiber-to-fiber transmission (Hi) 
Wavelength calibration (iv) Scientific image combination. 

We create a final data cube for each pointing by averaging over 
the final single-observing block (OB) flux-calibrated science cubes 
produced by the pipeline using the lEU tables, with the correspon¬ 
dence between each row (i.e. each spectrum) of the extracted spec¬ 
tra and the fiber’s position on the lEU head. The four lEU quadrants 
are treated independently up to the combination of all exposures 
into a single-pointing datacube. 

We subtract the median sky spectrum (extracted directly from 
the ‘sky frame’) from the ‘science frame’. Median combining was 
needed to ensure that any residual contamination from faint objects 
was removed. We are confident that this sky-subtraction procedure 
is accurate enough but, if prominent sky emission lines are still 
present in the spectra after sky subtraction, we mask them out dur¬ 
ing the the stellar kinematics recovery process. We also investigate 
possible effects of the differential atmospheric refraction by check¬ 
ing whether or not the photometric center of the galaxy changes 
its position as a function of wavelength in the datacube. We find 
that the position does not change by more than 2 pixels (~ 1.3 arc- 
sec) between 4000 and 6200 A. Each single-pointing datacube is 
astrometrically calibrated independently by matching the position 
of the photometric center or of a star (when the center is not visi¬ 
ble in that particular pointing) to an HST ACS image of the central 
region of NGC 4697. Subsequently, the different pointings are lin¬ 
early interpolated onto a common 3D grid uniformly sampled, after 
normalizing them for exposure times, using a simple grid interpola- 
tion[3to obtain the final mosaic datacube. As a further check of the 
robustness of this mosaic reconstruction, we re-check the WCS as¬ 
trometry on our final combined image, using as reference the same 
HST/ACS image. 

The results of this comparison are visible in Eigure which 
shows the final reconstructed mosaic image, obtained by averaging 
the datacube along the wavelength direction, as well as the ACS 
image. Twenty blue isophotes, derived from the HST image, are 
shown on both images to check the accuracy of our astrometrical 
calibration. Although most of the point sources visible in the HST 
image are too faint to be detected in our VIMOS observations, the 
positions of the center and of a star (ENE) coincide between the two 
instruments. On the right side of Eigurewe also plot a ID spec¬ 
trum extracted at the photometric center of the galaxy and another 
ID spectrum extracted at ~ 40”, in a position which is highlighted 
in the reconstructed EOV by a small yellow star. 


^ N. Devillard, Y. Jung, R. Paisa, C. Izzo, P. Ballester, C. Sabet, K. Banse, 
M. Kiesgen, L. Lundin, A. Modigliani, D. J. McKay, N. P. F. McKay, J. 
Moeller-Larsen, L. de Bilbao 

^ We interpolate, shift and rotate each single-pointing datacube by match¬ 
ing the position of the center and the rotational angle into this datacube- 
container. 
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Figure 1. Left-upper panel: NGC 4697 final reconstructed field-of-view generated with the mosaic of all the on-target pointings, summing the fluxes in the 
datacube over the whole wavelength range. Each dotted square represents a single pointing of 27" x 27". Green cubes are data from the first run (P79), 
whereas orange cubes are the two central pointings from the second one (P85). Left-bottom panel: HST ACS image of the central region of NGC 4697 showed 
at the same spatial scale as the reconstructed VIMOS image. Blue isophotes (derived from the HST image) are shown on both images to check the validity and 
precision of the astrometrical calibration (see text for more details). Right panels: ID galaxy spectra extracted in the central spatial bin (top panel) and in a 
spatial bin at ~ 40” from the galaxy center (middle panel), and ID sky spectrum (bottom panel) extracted from the central bin of the skycube. The spaxel from 
which the middle-panel ID spectrum has been extracted is highlighted on the reconstructed FOV by a yellow star. The spectra are plotted versus restframe 
wavelength and several spectral features of the galaxy such as G4300, H/3, Mg 6 and Fe lines are marked. Sky emission lines have been not removed from the 
spectra but have been properly masked out during the fit procedures. 


3 KINEMATICS MEASUREMENTS 


To recover the line-of-sight velocity distribution (LOSVD) first, we 
apply a spatial binning method in order to increase the signal-to- 
noise ratio (S/N) over the single-fiber spatially resolved spectra of 
the final datacube (Section 3.1). Then we measure the LOSVD for 
the galaxy using a maximum penalized likelihood approach (Sec¬ 
tion 3.2). In this way, we obtain spatially resolved 2D profiles 
of rotation velocity, V, and velocity dispersion, cr, as well as the 
higher order Gauss-Hermite moments hs and /i 4 , which quantify 
the asymmetric and symmetric departures of the LOSVD from a 
pure Gaussian (related to the skewness and kurtosis respectively). 
In the following sections we describe in more detail the spatial bin¬ 
ning procedure and the LOSVD recovery. 


3.1 Voronoi 2D spatial binning 

We re-bin the data with a fixed final S/N by locally averaging neigh¬ 
boring fibers (using the Voronoi Adaptive technique IDL code of 
[Cappellari & Copm||2003| l. We use an adaptive binning scheme 
where the size of the bin is adapted to the local S/N. Bigger spatial 
bins are created in the low-S/N regions, with higher resolution and 
smaller bins retained in the high-S/N regions (towards the galaxy 
center]^ 

We tested several values for the final Voronoi fixed S/N, 
choosing a final threshold of S/N = 70. In this way we are able to 
obtain a high resolution spectrum for each spaxel, which is funda¬ 
mental for performing a proper stellar population analysis based on 
line-index measurements, while also preserving spatial information 

Spaxels are arranged in a regular grid, so the size is a well defined quan¬ 
tity. 






















VIMOS mosaic of the ETG NGC 4697 5 


that is crucial for recovering detailed spatially resolved kinematics 
information up to ~ 0.7 effective radius 66 ”). 

3.2 Line-Of-Sight Velocity Distribution (LOSVD) fitting 

We measure rotation velocity, velocity dispersion and higher-order 
moments (V, cr, /i 3 , Ka) for each Voronoi-rebinned spaxel spec¬ 
trum using the pPXF code of |Cappellari & Emsellem| ( [2004) . This 
software, working in pixel space, finds the combination of stel¬ 
lar templates which, convolved with an appropriate line-of-sight- 
velocity distribution, best reproduces the galaxy spectrum. The 
best-fitting parameters of the LOSVD are determined by minimiz¬ 
ing which measures the agreement between the model and 
the observed galaxy spectrum, over the set of N good pixels. We 
use the MILES Library of Stellar Spectra ( |Sanchez-Blazquez et al.| 
|2006[[FSc6n-Barroso et al.|2011| : 985 well-calibrated stars cover¬ 
ing the region from 3525-7500 A at a spectral resolution of 2.54 
A (FWHM), cr ~ 57 kms“^ at the central wavelength, obtained 
at the 2.5m Isaac Newton Telescope (INT). The pPXF code also 
allows the user to mask out noisy regions of the galaxy spectrum. 
We mask out sky emission lines, focusing on absorption lines from 
4300 — 5400 A (including in this way H/3, O lines. Mg 6 dou¬ 
blet and Fe). Uncertainties on individual points are obtained via a 
Monte Carlo simulation. 

We follow the approach described in |Cappellari & Emsellem| 
( |2004| ), to recover rotation, velocity dispersion and the higher-order 
terms of the Gauss-Hermite series simultaneously. The method 
makes use of a penalization factor to the solution, which biases the 
results to suppress the noise. Obviously, one would like the penalty 
to leave the LOSVD virtually unaffected. Thus, we first perform the 
kinematics fit without including any penalty (i.e. setting the key¬ 
word BIAS=0 in the pPXF code) to make sure that our results have 
not been biased and also to properly estimate the uncertainties on 
the velocity distribution. This produces a noisy solution (especially 
for cr) but allows us to set a range of values for /13 and /14 which 
can be used to estimate the right value for the penalization factor. 

In particular, this has been done with the ‘ppxf simulation ’ 
procedure (included in the publicly available code). The routine 
performs a Monte Carlo simulation of the spectra, which repeats the 
full fitting process for a large number of different realizations of the 
data obtained by randomly adding noise to the original spectra with 
a resulting S/N equal to the chosen threshold value (S/N = 70 in our 
case), and using as input for [hs , Ha] the maximum representative 
values measured in the non-penalized pPXF fit of the previous step 
([—0.3,0.3]). At this point, we repeat the fitting procedure while 
using as penalty the largest value such that the mean difference 
between the output and the input [hs , /14] is well within the rms 
scatter of the simulated values (BIAS=1.5). 

Finally, we inspect by eye zones of the FOV where we expect 
to have contamination (e.g. near the brighter stars) and spectra for 
which the reduced of the fit was larger than 3. In particular, we 
identify ~ 30 “problematic” Voronoi cells, for which we go back to 
the single fiber spectra and re-create the Voronoi cell while exclud¬ 
ing contaminated or noisy spectra. Then we re-run the pPXF pro¬ 
cedure, in some cases also setting by hand the pixel regions where 
the fit is performed. 

We take into account all the main systematic uncertainties, 
mostly caused by template mismatch, and wavelength coverage. 
In order to minimize template-mismatch errors, we fit a broad set 
of stellar templates together with the kinematics, after convolving 
both to a common resolution. We also test the robustness of the 
kinematic results by changing the spectral region of the fit and we 


find results always consistent within 2a (not including systemat- 
ics). The uncertainties on the final inferred kinematics are estimated 
by adding in quadrature the formal uncertainty given by the Monte 
Carlo simulation and the scatter in the results for different templates 
and spectral regions. Typically, we find uncertainties of the order of 
5% for rotation and velocity dispersion and at most of the order of 
10% for the higher-order LOSVD moments. For the central bin, we 
obtain a velocity dispersion of cr* = 161 zb 8 kms“^ perfectly 
consistent with previously published results (i.e. |Binney, Davies"^ 
|Illingworth|1990[[de Lorenzi et al.|2008] l. 


3.3 Two dimensional kinematics maps 

The reconstructed mosaic (~ 80” x 60” effective area), is among 
the largest 2D maps ever measured for an elliptical galaxy and pro¬ 
vides a powerful tool to full characterize the internal kinematics of 
the system and the spatially-resolved properties of its underlying 
stellar population, which are necessary ingredients to understand 
the dynamics and its connection with formation mechanisms. Here 
below, we characterize the 2D behavior of all the relevant kinemat- 
ical quantities, while in the next sections we discuss the properties 
of the stellar populations and try to establish their connection with 
the kinematical structure. The full dynamical analysis as well as 
the derivation of the mass and orbital profile of the galaxy will be 
presented in a forthcoming paper. 

Figure shows 2D maps of rotation velocity, velocity disper¬ 
sion and the Gauss-Hermite moments hs and Ka of the LOSVD. 
On the left column, these quantities are shown as derived from the 
Voronoi-binning technique discussed in Section [3T] Here the small 
scale features are visible, mainly in the center, while the large scale 
structure at larger radii looks noisy. In order to make more evident 
the large scale pattern of the kinematical quantities, we show on the 
right column the 2D quantities smoothed by averaging the Voronoi 
values over a grid of 15 x 15 macro bins and then interpolated. The 
final average maps have been masked with an ellipse mimicking 
the orientation and the axis ratio of the outermost isophote enclos¬ 
ing the VIMOS mosaic. 

The galaxy shows a rotation velocity (upper panel) and hs 
(third panel from above) patterns typical of a rotation-supported 
structure along the major axis of the system. In particular, the hs 
indicates the presence of a disk-like structure elongated along the 
major axis. 

The velocity dispersion map (second panel), on the other hand, 
is rather fiat over the full area covered by our observations, with 
minimal, statistically insignificant fiuctuations. 

Finally the Ha parameter map (bottom panel) shows a clear 
elongated area along the line of nodes with positive values (Ha ~ 
0 . 01 ), which tend to become negative in the regions off the plane, 
which are not common for disk-like structure. The presence of neg¬ 
ative Ha off the main rotation axis in the regions where the rotation 
field is almost zero indicates tangential orbits, possibly belong¬ 
ing to the bulge component that starts to dominate outward (see 
Sec |4.1| ). The presence of two components along the line of sight 
might be responsible for the positive Ha parameter aligned along 
the galaxy major axis by producing a overall peaked distribution 
due to the superposition of a cold disk and hot bulge. 

IFS provides a complete and high resolution 2D map of the 
LOSVD, allowing the definition and estimation of a number of im¬ 
portant parameters that can be used to better understand the detailed 
internal structure of the galaxy, a crucial ingredient for investigating 
possible formation scenarios. An interesting quantity that is possi- 
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Figure 2. NGC 4697 2D line-of-sight velocity distribution maps in a 80” x 
80” box. Panels on the left column show the un-smoothed Voronoi-rebinned 
version, as derived in Sec. |3.1| Here the small scale features are visible. 
Panels on the right column show the 2D maps smoothed by averaging the 
Voronoi values over a grid of 15x15 macro bins and then masked with 
an ellipse mimicking the orientation and the axis ratio of the outermost 
isophote enclosing the VIMOS mosaic. Here the large scale pattern of the 
kinematical quantities is more appreciable. 


ble to derive in 2D is the ratio of the ordered and random motion in 
a galaxy (V/a). 

In Figure]^ we present the first V/a 2T) map of NGC 4697, 
up to ~ 0.7i7eff,which gives a more quantitative view of the dy¬ 
namical state of the galaxy and further confirms the presence of 
the galaxy sub-components. One expects that the regions which 
are bulge-dominated have smaller V/a values, while the rotation- 
supported regions where the disk is dominant have larger V/a val¬ 
ues. This is indeed what we observe in the 2D map and it is further 
confirmed in the next Section, where we compare our /is- and /14- 
V/a correlations with those inferred from hydrodynamical simula¬ 
tions. 

We note that none of the previous 3D measurements of galax¬ 
ies have performed the same kinematical diagnostic as the one 
above. However, more standard rotation and dispersion 2D maps 
of the LOSVD of NGC 4697, like the ones derived from the 
ATLASsurvey, are fully consistent with the ones presented here, 
although with a lower spatial resolution and over a smaller FOV 



AX[arcsec] 


Figure 3. Top panel: V/a 2D Voronoi map, quantifying the ratio of the 
ordered and random motion of NGC 4697. 


(ATLAScovers 40 x 60 arcsec, Krajnovic et al.|201l| . We do 
not directly compare our 2D maps to those obtained by ATLAS 
because of the different spatial and spectral resolution but in Sec¬ 
tion 3.5 we extract ID profiles along the major axis and compare 
our kinematics profiles against the ATLAS^^ones. 

Finally, extended 2D kinematical maps up to 3.5 i^eff of this 
galaxy were recently presented as part of the SLUGGS Survey 
( [Arnold et al.|2014] ). Overall there is a fair agreement on the fea¬ 
tures observed in the data, both in the velocity field and in the 
higher-order moments of the LOSVD. In fact, based on their more 
spatially-extended kinematic map, the SLUGGS team concluded 
that the fast rotation is due to an embedded disk with a rotation am¬ 
plitude clearly declining with radius. The presence of a disk is also 
found from our detailed analysis of the /13— and h^-V/a diagrams. 


3.4 Comparison with simulations 

hs-V/a and h^-V/a correlations have also been predicted in sim¬ 
ulations of disk merging ( jNaab, Jesseit & Burkert|200^ H-i-09 and 
N-i-14). For instance, N-i-14 linked the present day central kinematic 
properties of galaxies to their individual cosmological formation 
histories and showed that galaxy mergers (major and minor) have 
a significant influence on the rotation properties, resulting in both 
a spin-down and a spin-up of the merger remnants, and in a clear 
anti-correlation between /13 and V/a. 

N-I-14 have defined six classes (named with letters from A to 
F) with different evolutionary paths and histories of mass assem¬ 
bly that can be directly linked to their 2D kinematic maps and the 
present day shapes. The galaxies groupings are based on the 2D 
kinematic maps and on the so-called ‘spin’ parameter (Ar), a proxy 
for the protected angular momentum per unit mass (jEmsellem et al.j 
|2007l|20n] ). 

The Ar parameter has been used to show that, within the ef¬ 
fective radius, a dichotomy exists between ‘fast-’ (Ar> 0.1) and 
‘slow-’ (Ar< 0.1) rotators ( jEmsellem et al.|[2007] ). The second 
class of galaxies shows little or no rotation, significant misalign¬ 
ment between the photometric and the kinematic axes, and contains 
kinematically decoupled components, whereas the first class ex¬ 
hibits regular stellar velocity fields, consistent with disk-like rota- 
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tion and sometimes bars. This version of the ‘spin’ parameter takes 
the form: 




( 1 ) 


where Fi, Ri, Vi and ai are the flux, circular radius, velocity and 
velocity dispersion of the i-th spatial bin, with the sum running over 
the Np bins. 

A revised and more accurate separation of these two classes 
of galaxies has been proposed by |Emsellem et al.| ( |2011| ) to include 
the associated apparent ellipticity: a higher value of the specific 
stellar angular momentum is expected for galaxies which are more 
fiattened or closer to edge-on if these are all intrinsically fast rota- 

tor|3 

The threshold for Ar was fixed to be proportional to square- 
root of £ with a scaling parameter that depends on the aperture (see 
Appendix B in |Emsellem et al!]|2011| for further details). For in¬ 
stance, at i^eff /2 


= (0.265 ± 0.01) X ^/(ee///2) (2) 

Following this classification, NGC 4697 was previously clas¬ 
sified as “fast-rotator”. This means that, given an ellipticity of £ = 
0.36 ( |Peletier et al.|1990] l, one expects for this parameter a value 
of Ar> 0.16 (for r < Res/2). We find Ar ^^,/2 = 0.35 ± 0.05, 
confirming the classification of this galaxy as a ‘fast-rotator’ and 
in perfect agreement with the ATLAS^^result (Ar ^^^/2 = 0.322, 
[Emsellem et al.|2011| ). 

Fast-rotators belong to class A, B and D in the N-i-14 classi¬ 
fications. However only the first two groups of galaxies will result 
in hs and distributions similar to that we observe in NGC 4697. 
Both groups (A and B) of galaxies are rotationally supported, had 
a significant amount of central in-situ, dissipative, star formation 
but, while class A have late assembly histories dominated by minor 
mergers, class B galaxies have gone through a late gas-rich major 
merger. These two formation scenarios will produce the same hs- 
V/a anti-correlation that we observe in our maps. 

One way to distinguish between the two classes is to look at 
the Ar profiles: class A galaxies show peaked Ar profiles within 
the effective radius, originating from the fast rotating central stellar 
disk, whereas class B galaxies have Ar profiles constantly rising 
beyond the effective radius. 

In Figure we plot the Ar profile as a function of radius 
calculated from circular apertures. The profile has a strong peak 
around 0.2i^eff and then it rises until , which suggests that 
the mass assembly of NGC 4697 has involved significant in-situ 
star formation, that NGC 4697 has experienced several gas-rich mi¬ 
nor mergers (like simulated galaxies of class A) and possibly a late 
gas-rich major merger leaving a previously rapidly rotating sys¬ 
tem unchanged (like simulated galaxies of class B). However given 
the poor resolution of the simulations this interpretation should be 
taken with caution (Fig.5 in N-i-15 shows the Ar profiles of the 
simulated galaxies in the six classes and only one galaxy of class 
A shows a peak, which is at ~ 0.95i^eff- On the other hand we 
observe that in Fig.5 of |Emsellem et al.|2011| showing Ar for the 


^ A galaxy with a relatively low value of Ar, e.g., of 0.2, may be consistent 
with a simple spheroidal axisymmetric system viewed at a high inclination 
(near face-on), but this is true only if its ellipticity is correspondingly low 
(e.g., £ ^ 0.3). A large ellipticity would instead imply a more extreme 
object (in terms of orbital structure or anisotropy). 



Figure 4. The Ar parameter calculated using the formula of |Emsellem| 
|et al.|j2007) as a function of radius. 


full ATLASsample, several fast-rotators show a profile simi¬ 
lar to what we derive for NGC 4697). Moreover, |Coccato et al.| 
\2009\ derived the Ar profile for NGC 4697 out to ~ 3.5i^eff 
based on planetary nebulae (PNe) and showed that it flattens out 
at ~ which supports the idea that NGC 4697 should be 

compared with simulated galaxies of class A. We also note that 
[Arnold etliLlpOMI ) concluded that galaxies like NGC 4697 formed 
through rapid, dissipative processes, which formed the central re¬ 
gions, followed by a prolonged period of dry minor mergers. We 
cannot exclude the scenario of gas-rich late major mergers, but in¬ 
dependent lines of evidence point toward a scenario where the mass 
assembly of NGC 4697 is dominated by minor mergers and gradual 
dissipation. 

In Figure|^we show the hs- and h^-Vja correlations derived 
from the 2D kinematic maps (left panels), and we compare them to 
luminosity contours derived by H-i-09 for a 1:1 disk-disk merger 
with 30% gas fraction (middle panels) and from the cosmological 
hydrodynamical simulations of N-i-14 (class A and class B, right 
panels). A good agreement, at least on a qualitative level, is found 
between data and simulations, although we caution that H-i-09 only 
considered the case of equal-mass mergers. 

VIMOS data-points are color coded in accord with their value 
of |Cos[^]|, calculated with respect to the galaxy minor kine- 
matical axis. We find a clear correlation between the |Cos[^]| 
values and the regions where the disk component is elongated. 
Points oriented near the major axis (i.e., |Cos[^]| ~ 0) have also 
higher absolute values of /i 3 and therefore are associated with the 
rotation-supported structure. Similarly, points with higher values of 
I Cos [0] I have lower /i 3 and roughly correspond to the bulge region. 
The /i 4 does not show a similar (anti-)correlation with V/a. 

The h‘i-V!a anti-correlation generically tracse cold disks. It 
is a feature that naturally emerges in disk-disk gas-rich mergers, 
and tends to produce large asymmetric deviations from Gaussian- 
ity because of the presence of sub-components that retain a mem¬ 
ory of the disk-star streaming, and because of gas that settles down 
through dissipation to form a thin embedded disk. This hs-Vja 
correlation is generally negative when the dominant orbit popu¬ 
lation tends to stream (e.g. z-tube orbits, see also [Rottgers, Naabj 
|& Oser|[2014| ), and positive when the rotation is caused by sub¬ 
dominant streaming population (e.g. x-tube orbits) superimposed 
on the tail of the dominant streaming orbits (e.g z-tubes which still 
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Figure 5. Left panels:The hs-Vja (top) and the Ka-V j(J (bottom) distributions from the VIMOS kinematic maps. Points are color-coded according to the 
azimuthal deviation form the galaxy minor axis. Middle panels: Comparison with predictions from 1:1 disk-disk merger simulations from H-i-09 with a 30 
percent gas fractions. The shading is proportional to the logarithm of the luminosity, whose scale is plotted in the bar on the top and the red contours contain 
68% and 95% of the luminosity. Right panels: Comparison with predictions of N-i-14 simulations for galaxies of class A and B. 


contribute some rotation along the minor axis; see also the discus¬ 
sion in H-i-09). The same features are not found in second genera¬ 
tions of major (re-)mergers (mergers at lower redshift, where less 
cold gas is present and with elliptical-morphology progenitors), in 
which case the features are more mixed and look more concentrated 
around the center of the diagram. 

The hA-VI a correlation of NGC 4697 in Figure shows a 
V-shaped distribution with a double peak at \V/a\ ~ 0.5. This 
distribution is very similar to the predictions from simulations of 
N-i-14 (right, bottom panel) and it also agree with H-i-09 (middle, 
bottom panel), although there is a significant offset between the 
data and the H-i-09 predictions. The peaked population could cor¬ 
respond to a late-formed disk stars (note the correspondence to the 
extreme hs values and | Cos [0] \ ~ 0 indicating the alignment along 
the major axis) which tend to produce a positive Ha in superpo¬ 
sition with the more isotropic component (differences in ages of 
the two components are discussed in Sec.[^. In the simulations of 
H-I-09 the re-mergers would look more concentrated around the ori¬ 
gin and would tend to have LOSVDs with positive Ha, similarly to 
expectations from violent relaxation, which generally results in a 
radially-biased distribution function ( [Binney & Tremaine|2008| ). 

To conclude, the high-order moments of the LOSVD provide 
kinematical evidence for a two-component system: the largest ab¬ 
solute values of hs and Ha indicate an area of the galaxy with 


the dominance of x-tube over z-tube orbits, which are produced in 
collisionless mergers and trace the presence of a (later) streaming 
structure, these are also the regions around |Cos[^]| ~ 0 (minor 
axis). Values with /is ~ 0 and /14 ~ 0 indicate instead the presence 
of a non rotating subsystem that dominates in the regions where 
the values of \Cos[0] \ are higher. The hs-Vjcj and the peaked Ar 
profile suggest that gas-rich minor mergers dominate the mass as¬ 
sembly of NGC 4697, although we cannot completely exclude the 
possibility of a major merger event that happened at redshift \. 

We make use of this evidence in Section [TT] where, based on 
the photometry of NGC 4697, we separate the disk from the bulge 
and study the ages and metal abundances of the stellar populations 
of these two components. 

The details of the gas accretion leading to the formation of 
the central disk may be somewhat different than the one investi¬ 
gated in the simulation of disk merging and it might be useful to 
have LOSVD predictions from simulations of cold accretions as an 
alternative high dissipative process forming the disk. It seems evi¬ 
dent, though, from H-i-09 and N-i-14 that dry mergers is not able to 
produce the kinematical substructures we see in the /i3- hL-Vja 
diagrams. In fact, in a merger without a dissipative gas component, 
the remnant galaxies are axisymmetric and do not show a popula¬ 
tion of stars with high angular momentum tube orbits nor the re¬ 
growth of a significant disk component (e.g. class D in N-i-14). 
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Figure 6. NGC 4697 ID kinematic profiles along the major axis. Black 
points are da ta from this work, blue points are long-slit data from |de Lorenzi| 
|et al.||2008) and green points are IFU data from the ATLASsurvey, ex¬ 
tracted along the major axis with the same procedure used to extract our 
points. A good agreement is found between the different studies. Uncertain¬ 
ties on the rotation and velocity dispersion are of the order or 10 km s“^, 
while on the hs and /i 4 parameter are of the order of 10 percent. 


3.5 One dimensional spatially resolved kinematics profiles 

To compare our results against published results based on long-slit 
data and IFS with different FOV and resolution, we extract from the 
2D kinematics maps the information along the major axis. Specif¬ 
ically, we extract and sum all the spectra positioned in a rectangle 
with width of 1.34" (2 pixels) oriented along the major axis of 
the galaxy and centered on our kinematic center (to be as close as 
possible to the configuration of slit-based studies). We perform the 
same extraction for the IFU ATLAS^^observations (publicly avail¬ 
able on the Survey Webpage) in order to perform a simple com¬ 
parison of the two sets of data with different spatial and spectral 
sampling. In Figurewe compare our profiles (black points) with 
the ATLAS^^ones (green points) and with the ones published in 
de Lorenzi et al. (2008, hereafter DL08), derived from high S/N 
integrated absorption line spectra obtained by [Mendez et al.| ( [2005| ) 
with FORS2 at the VLT, a slit width of 1 arcsec and seeing 1-1.5 
arcsec (blue points). For the major axis inclination, we assume the 
value PA = 63 degrees from jBinney, Davies & Ilhngworth|(r990j ), 
slightly different from the one assumed by DL08 but providing the 
best symmetry and maximum rotation for our extracted ID pro¬ 
files. We do not worry about this small PA discrepancy that could 
be caused by the fact that the slit in DL08 was not exactly on-center 
or could be related to the offset we found between kinematic and 
photometric center, most probably caused by a small effect of at¬ 


mospheric refraction. For the ATLAS^^ data we use the same in¬ 
clination of PA = 63 degrees and extract the major axis profiles 
with the same procedure applied to our VIMOS data. We found an 
overall good agreement between the three ID kinematical profiles 
which nicely demonstrates the robustness of our mosaic approach 
and the validity of our results. 


4 2D STELLAR POPULATION ANALYSIS 

In this section we study the stellar population of the disk and of 
the bulge of NGC 4697 separately, to put constraints on the age, 
[a/Fe], metallicity and IMF normalization of the stellar population 
of the different morphological components of the galaxy. 

To constrain the stellar population parameters, we focus on 
line-index measurements rather than full spectral fitting, which 
avoids issues with spectral calibration when comparing to obser¬ 
vations that might have been poorly calibrated and we assume a 
simple stellar population (SSP), rather than an extended star forma¬ 
tion history, which is still a strong assumption and will be properly 
addressed in future work. 

We measure equivalent widths (EWs) of several stellar absorp¬ 
tion features in each of the Voronoi-rebinned ID spectra, after con¬ 
volving them to a common resolution of 350 kms“^ to correct 
for kinematic broadening We compared these EWs to those mea¬ 
sured from an extended version of the SSP models of Conroy & van 
Dokkum (2012a, hereafter CvD12) with the same common resolu¬ 
tion, obtained in Spiniello, Trager & Koopmans (2014b). 

CvD12 optimized models over the wavelength interval 
0.35/xm < A < 2.4/im at a resolving power of R ~ 2000, 
specifically for the purpose of studying old, metal-rich stellar pop¬ 
ulations. The models are built combining two empirical stellar li- 
braries (MILES, [Sanchez-Blazquez et'^|2006[ and IRTE, jCush-j 
jing, Rayner & Vacca||2005| l and three sets of isochrones to mea¬ 
sure the IME slope down to ~ 0.1 Mq. The CvD12 models ex¬ 
plore variations in age in the range 3-13.5 Gyr, a—enhancement 
of 0-0.4 dex, and four different single-slope IMEs: a bottom-light 
single-slope power law with x — 1.^ (where x is the IME slope, 
using dN/dm oc m~^) that produces a stellar mass-to-light ra¬ 
tio comparable to the one obtained from a standard Chabrier IME 
( |Chabrier|2003j ), a Salpeter (1955) IME with a slope of x = 2.35, 
and two bottom-heavy (dwarf-rich) IMEs with slopes of x = 3.0 
and X = 3.5. 

The CvD12 models also allow for variations in the [a/Ee] ra¬ 
tio as well as variations in the abundance pattern of 11 different 
single elements. This is enough to decouple IME variations from 
abundance variations. However, all of the CvD12 models use so¬ 
lar metallicity isochrones, even when synthesizing with different 
abundance patterns or different [a/Ee]. The abundance variations 
of single elements are implemented at fixed [Ee/H], which implies 
that the total metallicity Z cannot be easily controlled. This is re¬ 
strictive, if one aims to disentangle elemental enhancements, metal- 
licity changes and IME variations, especially in the case of massive 
ETGs, which have star formation histories different from the solar 

® This final resolution has been chosen for practical reasons: the version of 
the simple stellar population models used in this paper was created and used 
in S+14 with a resolution of 350 km . To save computer power, we did 
not perform again the EW measurements and simply smoothed our data to 
the resolution of the models. We are confident that this did not cause any 
loss of information, and we check (by eye-inspection) that contamination 
from sky line residuals do not affect the spectral features. 
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neighborhood (e.g., |Peterson||197^ |Peletier|1989t |Worthey|l^ 


[Trager et al.|2000b||Arrigoni et al.|2010ir 

To resolve this issue, we follow the approach presented in 
[Spiniello, Trager & Koopmans] ( |2015] ), hereafter S+15 of extend¬ 
ing the parameter space of the SSP models using response func¬ 
tions. We construct a metallicity-response function for each given 
age and IMF slope taki ng the ratio between t wo spectra of the 
MIUSCAT SSP models ( Vazdekis et al.|201^ with the same age 
and IMF slope but different total Z. We then multiply this re¬ 
sponse function with the spectrum of a CvD12 model with fixed 
age and IMF to build a new model (SSP) that extrapolates the lat¬ 
ter model to a new part of parameter space (i.e. covering super¬ 
solar Z). These modified SSP models combine the fiexibility of 
both the MIUSCAT and CvD12 models to predict spectra in a part 
of parameter space that neither of them reaches separately (MIUS¬ 
CAT models use different total metallicities but do not allow one 
to change the abundances which are fixed to solar). As pointed out 
in S-i-14 and S-i-15, selecting IMF-dependent features that are age- 
and metallicity-independent and combining them with indices that 
depend mainly on age or mainly on element abundance is impor¬ 
tant for breaking the age-metallicity-IMF degeneracy when using 
SSP models to infer the stellar populations from unresolved galax¬ 
ies spectra. 

We measure some of the classical blue Lick indices (e.g. H^, 
Mg6, Fe5270, Fe5335, [MgFe]: |Gonzalez||1993[ |Worthey|[T9^ 
[Trager et al.pOOOat that have been extensively and successfully 
used in the literature to constrain the age and metallicity of galax¬ 
ies, together with some of the bluer IMF-sensitive features defined 
and used in S-i-14 and S-i-15. Unfortunately this paper makes use of 
data taken with the HR-blue grism, and therefore only covers the 
bluer IMF-features of S-f-14 (specifically bTiO, aTiO, TiOl) while 
in the future we plan to obtain new data using the HR-red grism 
which will cover also Ti02, CaHl, Ti03 and Ti04, allowing a more 
detailed analysis and a more secure inference on the low-mass end 
of the IMF slope. 

Before measuring EWs, we convolved galaxy and model spec¬ 
tra to an effective velocity dispersion of a = 350 km to correct 
for kinematic broadening. Indices in both the galaxy and the model 
spectra have been measured with the same definitions and meth- 
odfl We do not place our indices on the zero-point system of the 
Lick indices and quote them as index strengths in units of angstrom 
(if atomic) or magnitudes (if molecular, like the TiO lines). 

The great advantage of integral field spectroscopy is the si¬ 
multaneous availability of spatial and spectral information. This 
allows us to generate 2D spatial maps of spectral features of in¬ 
terest, like the one in Figure]^ showing the [MgFeQindex that has 
been widely used to constrain the metallicity. The 2D map suggests 
different stellar populations in the bulge and in the disk: specif¬ 
ically, the [MgFe] EWs are systematically higher in the Voronoi 
cells within the disk. Therefore, the 2D element maps indicate that 
the stars in the disk are more metal-rich than the stars in the bulge. 
However, we caution the reader that it is not possible to directly in¬ 
fer the total metallicity by looking at a single-element map, because 
of the degeneracies present between the stellar population param¬ 
eters. Eor instance, in this case the [MgEe] values could be higher 


^ The MIUSCAT SSP has instead a fixed total metallicity and covers six 
metallicity bins Z = 0.0004, 0.001, 0.004, 0.008, 0.019 and 0.03, where 

0.019 represents the solar value. __ 

® We use the code SPINDEX from Trager, Faber & Dressier |2008 i. 


[MgFe] = v/(Fe5270 + Fe5335)/2 x Mgb, {ConUe^ 
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Figure 7. NGC 4697 2D maps of the [MgFe] stellar absorption feature mea¬ 
sured from the spectra at a final resolution of a = 350 km The map 
clearly shows that the stellar population of the disk is different from the 
stellar population of the bulge component of the galaxy. Areas within the 
disk are infact more metal-rich (higher [MgFe] EWs). 


due to an increase of age and/or could increase of total metallic¬ 
ity. This is an example of the so-called age-metallicity degeneracy 
which can be broken by using more than one indicator. This moti¬ 
vates us to decompose the bulge and disk components and to extract 
a ID spectrum of each of them in order to study their stellar popu¬ 
lations separately through a more quantitative, proper SSP analysis. 


4.1 Disk-Bulge decomposition 

We perform a two-component fit to the HST V-band image of 
NGC 4697 obtained by |Tal et aL| ( |2009| ), driven by the correlations 
in the hs and parameters versus V/a that clearly indicate the 
presence of a disk. We use the latest version of GALFIT ( |Peng| 
|et al.|2010] | to fit two ellipsoid models (‘sersic’ for the bulge and 
‘expdisk’ for the disk) to the light profiles of the galaxy image. We 
obtain a similar position angle (PA ~ 64 deg, confirming previ¬ 
ous findings) but very different axis-ratios for the two components: 
{b/a) disk = 0.297, {b/a)buige = 0.7512. We use those to spec¬ 
troscopically separate the bulge from the disk. We note that a more 
rigorous spectroscopic separation, following the approach of |Coc-| 
|cato et airi ( |2011| ) and [Coccato, lodice & Arnaboldi| ( |2014^ , will be 
performed in future studies to fully separate these two components. 

First, we extract and sum together all the spectra in the dat- 
acube with r <i?eff/8. In this central region we expect bulge and 
disk to both contribute substantially. Second, we select and sum the 
spectra belonging to an elliptical anulus centered on the kinematics 
center with the axis-ratio and PA found by GALFIT for the disk: 
{b/a)disk = 0.297, PA = 65.47deg and with r >i^eff/8 (to ex¬ 
clude the central region). Finally, we sum up the remaining spectra 
where the contribution of the bulge is dominant. The resulting ID 
spectra have very high S/N (S/N = [450,320, 280] for center, disk 
and bulge respectively). 

In Figurewe show the extracted ID spectra of the disk and 
bulge, as well as of the mixed (disk+bulge) population of the cen¬ 
tral region of NGC 4697 for comparison. Shaded yellow bands 
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Figure 8. ID spectra extracted from the bulge (blue), disk (red) and the central /§ region (10 arcsec, black). Yellow boxes show the band-passes of the 
line indices used to constrain the SP paramenters, while blue and red bands represent the blue and red continuum side-bands respectively. 


show some of the line indices used in this work, along with contin¬ 
uum side-bands (blue and red colored bands). Some small differ¬ 
ences are visible already from the spectra: the disk spectrum shows 
stronger Mg and Fe lines, which confirms what have already seen 
from the [MgFe] 2D map and gives hints of a more metal-rich pop¬ 
ulation. The NaD absorption line, the narrow strong line around 
5900 Ais stronger in the bulge than in the disk, consistent with the 
hypothesis that bulges could be Na-enhanced (e.g. in the bulge of 
the Milky Way, on average [Na/Fe] ~ 0.2, |Fulbright, McWilliam ^ 
|Rich|2007] ). However, we do not use this particular line in the fol¬ 
lowing analysis because its blue band is heavily contaminated from 
emission sky-lines and because of the recent finding presented in 
S-i-15 showing that the NaD lines are model dependent and there¬ 
fore must be interpreted with caution. The aTiO index has also been 
excluded from the analysis because it is contaminated from emis¬ 
sion sky-lines. 

In Figure we show index-index plots of the selected optical 
indicators to provide qualitative constraints on the stellar popula¬ 
tion parameters. In panels a) and b) the IMF dependence is minimal 
and the age and abundances of the galaxy can be inferred. The IMF 
normalization of the galaxy can be instead inferred from panels c) 
and d), using SSP models with the metallicity and age previously 
determined. 


4.2 Age and metallicity of the stellar population 

In the upper panel of Figure we show index-index diagrams of 
some of the classical Lick-indices. These diagrams are widely used 
to constrain age and metallicity of integrated stellar populations. 
The average indices extracted from the different regions are plot¬ 
ted on top of grids of CvD12 SSP models with varying ages (blue 
lines, [7,9,11,13.5] Gyr), [a/Fe] (red lines, [0.0, 0.2, 0.4] dex) and 
total metallicity (green lines, [-0.22, 0.0, -1-0.22]). All plotted mod¬ 


els have a Chabrier-like IMF, but for these indicators (H^, Mg6, 
Fe5270, Fe5335) the dependency on IMF is minimal. 

Panel a) of Figure shows the H/3-[MgFe] diagram. The H/3 
index is one of the most used age-indicators in the optical, and al¬ 
lows one to break the age-metallicity degeneracy, when combined 
with the commonly-used [MgFe] index. However, this diagram 
alone does not allow one to also break degeneracies with [a/Fe]. 
This is clear for instance in the case of the disk component. From 
panel a) it is inferred that the bulge has an old population with solar 
[a/Fe] and slightly sub-solar metallicity. The central region is also 
old but has super-solar abundances. For the disk the interpretation 
is somewhat more difficult. In fact, the disk spectrum matches a 
model with age = 11 zb 2 Gyr, slightly super-solar metallicity and 
solar [a/Fe] but also a model with older age, lower metallicity and 
a-enhanced. To discriminate between these two possibilities, the 
reader should refer to the results plotted in panel b). Panel b) shows 
Mgb versus < Fe >, which is the mean of Fe5270 and Fe5335. 
Here the disk matches only with a model with age = 11 zb 1 Gyr, 
solar metallicity and solar [a/Fe]. 

Index-index plots are a useful tool to give qualitative infer¬ 
ence on the stellar population parameters, but it is sometimes hard 
to break degeneracies and quantitatively constrain age, metal abun¬ 
dance and IMF slope at the same time. Therefore, in order to give 
more quantitative constraints on the stellar population parameters 
of the three components, we perform a minimization comparing 
the ID spectrum of each component with grids of linearly interpo¬ 
lated SSPs spanning a range of ages (between 7 and 13.5 Gyr with 
steps of 1 Gyr,), [a/Fe] (between 0.2 and -1-0.4 dex, with a step 
of 0.1 dex), and total metallicity ([M/H]=[-0.22,0.0,-i-0.22]) for a 
Chabrier IMF (we note that a moderate change in the IMF slope 
has only a negligible effect on these blue features). We find that the 
disk is slightly younger and more metal-rich than the bulge, con¬ 
sistent with a SSP model of ~ lO.Sl^'o solar total metallicity 
([M/H]=-0.03l°;f) and solar [a/Fe]([a/Fe] = +0.04l°;°2). The 
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Figure 9. Index-index plots of some of the classical stellar absorption features and of the bluer IMF-sensitive optical indicators defined in S+14 in the optical 
regime. Grids are CvD12 extended SSP models obtained in S+15 via the response function approach (see text for further details) . Blue lines show models 
with varying ages with values [7,9,11,13.5] Gyr, red lines are SSP models with varying [a/Fe]: [0.0,0.2,0.4] dex, green lines are SSP models with varying 
total metallicity ([M/H] = [—0.22, 0.0, +0.22]) and, finally, black lines are models with varying IMF slope from x = 1.8 (Milky Way-like) to x = 2.6 
(being Salpeter x = 2.35). Colored arrows in the corner of each panels show the direction of the variation of these three single parameters. Points with error 
bars are high S/N galaxy spectra extracted from the different regions by summing the Voronoi-rebinned spectra. Galaxies and models are convolved to a final 
common resolution of 350 km before measuring index strengths. In the plots on the upper row (panels a and b), IMF dependence is minimal (all models 
have x=l .8 here) and the age and abundances of the galaxy can be inferred. The IMF normalization of the galaxy can be instead inferred from panels c) and d), 
using SSP models with the metallicity and age constraints previously determined. Here, lines are CvD12 models with varying IMF slope in the range 1.8-2.6. 


bulge fits better with a 13.51 J '4 Gyr model with sub-solar metallic- 
ity([M/H]= -O.IT+J^) and solar [a/Fe]([a/Fe]= +0.03lo:oi)- 
The central region shows stellar population properties in between 
those of the disk and those of the bulge, as expected. It is as old 
as the bulge (12.9118 Gyr), shows a solar or slightly super solar 
[a/Fe] ([a/Fe]= +0.091 q qi) and it is even more metal enriched 
([M/H]=+ 0 . 231 q;q 5 ). These SSP results nicely confirm the qual¬ 
itative picture obtained via index-index plots in the upper row of 
Figure and are summarized in Tab. However, we note that 
since the ID spectra of bulge and disk are composite spectra where 


the contributions of the two components are not perfectly disen¬ 
tangled, the true stellar population differences are likely to be even 
greater. 

Finally we consider NGC 4697 as whole and calculate gradi¬ 
ents in its stellar populations. Following the approach of |Tortor^ 
|et al.| ( [^ 10 a| ), we define the stellar population gradient as the di¬ 
mensionless coefficient of the relation 

logX --logR/RM^x = (3) 
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Table 1. Best-fit stellar population models for the different components ob¬ 
tained via a ht with CvD SSP models with a x=1.8 IMF slope. 


Component 

Age (Gyr) 

[M/H] 

[o/Fe] 

Central 

Bulge 

Disk 

i2.9+i;| 

+o.23+^;^« 

-0.03lo;“ 

+0.09+^;^^ 

+0.03l°;“ 

+0.041°;“ 


where X represents the age or the total metallicity expressed in 
terms of [M/H]. We find a non negligible but very small age gradi¬ 
ent and a negative metallicity gradient in the sense that NGC 4697 
is more metal-rich and older in the center. These gradients are 
consistent with the measured gradients in old early-type galaxies 
( [Kuntschner et aL|2010||Tortora et al.|2010a| l and with remnants of 
gas-rich major mergers from simulations (e.g.|Hopkins et al.|2009| 
H-h09). 


4.3 The IMF normalization 

Using the metallicity and age constraints obtained from the Lick 
blue indices via the minimization, we focus now on two of the 
IMF-sensitive optical indicators that fall in the VIMOS wavelength 
range, the bTiO and the TiOl. We plot them against the [Mg/Fe] 
index in the lower panels of Figure The combination of the lat¬ 
ter indicator with IMF-sensitive features breaks the degeneracies 
between the [a/Fe], total metallicity and IMF slope variations. 

For each component, we use SSP models with varying IMF 
slope (x between 1.8 and 2.5) and [a/Fe], and with the other stellar 
population parameters as close as possible to the ones derived from 
the x^ fit to the blue indices. For the bulge, we use SSP models with 
13.5 Gyr and [M/H] = -0.22, for the disk we use a SSP models with 
11.0 Gyr and [M/H] = 0.0 (solar). Finally, for the mixed, central 
population we use SSP models with 13.5 Gyr and [M/H] =-i-0.22 
(super-solar). 

In Figure|^black lines show models with different IMF slopes 
and red lines show models with varying [a/Fe ]. No detectable vari¬ 
ation in the IMF normalization between bulge and disk is visible 
from the two panels. Both components of NGC 4697, as well as the 
central region, are consistent with a normalization of x = 2 . 0 ± 0 . 2 , 
slightly dwarf richer than a Chabrier IMF (x=1.8, still consistent 
within 1-sigma) but shallower than a Salpeter (x=2.35), which is 
perfectly consistent with the slope predicted from the IMF-cr,^ rela- 
tion (e.g.,|La Barbera et al.|2013[|Tortora, Romanowsky & Napoli-| 

tano|2()13[ S-fl4) for an ETG with = 160 kms“ . IMF nor¬ 

malization does not change with radius for this galaxy is consistent 
with the results found in |Martm-Navarro et al.| ( |2015] ). The lack of 
IMF spatial gradient is also consistent with the results on the low- 
mass system in |Martm-Navarro et~^ ( |2015| ) that shows an IMF 
normalization constant with radius and consistent with a Kroupa- 
like profile. 

4.4 Linking kinematics and stellar populations 

The main conclusion from the stellar population analysis is that the 
inner region (typically ~i?eff/ 8 ) looks coeval (or a bit older) and 
more metal-rich than the regions farther out, which suggests that 
this central region has been formed early via rapid collapse. How¬ 
ever, from the stellar population analysis alone, we cannot exclude 
minor or major merging of late-type disk galaxies, which have de¬ 
stroyed spiral arm features and left a younger disk component. It is 


therefore necessary to connect and correlate the stellar population 
results with those on the spatially resolved kinematics. 

The SAURON team ( [Kuntschner et al.|2006| jKrajnovic et al.j 
|2008| l have found that the fast-rotating component almost always 
features a higher metallicity compared to the galaxy as a whole. 
They have attributed this to the presence of a later-forming stel¬ 
lar population in a rotating structure with disk-like kinematics 
( [Kuntschner et al.|[2010| l. From the NGC 4697 2D maps as well 
as from the stellar population analysis presented in this study, we 
infer the same scenario. This is demonstrated inFigure which 
shows the correlation between the hs (top panels) and the (bot¬ 
tom panels) parameters and the FL/3 (left) and Mgb (middle) values 
obtained from the smoothed 2D maps. We run a Spearman test on 
these relations to better quantify the statistical dependence between 
the high-order moments of the LOS VD and FL/3 and Mg b respec¬ 
tively. The resulting Spearman rank correlation coefficients are p = 
0.36 (0.2) and p = 0.32 (0.36) for /i 4 (fe) and H^ and Mgb re¬ 
spectively. Given the big sample size (more than 100 points), these 
values of p all exceed the critical value {pcrit{n > 100) = 0.165) 
and therefore the relations are statistically significant, although not 
strong. 

The right panels of Figure [Toj shows the H/3 -h^ (top panel) 
and -/i 4 (bottom panel) diagrams with the points color-coded ac¬ 
cording to the Mg5 values. Although the discussion must be kept 
on a qualitative level because of the well-known age-metallicity 
degeneracy ( [Worthey|1994[ ), we can simplify the situation and say 
that larger H/3 values indicate younger ages and larger Mgb values 
indicate higher total metallicity. Thus we can conclude that points 
with /i 4 > 0 , tracing the later-formed disk, are slightly younger 
and more metal rich, while points with /i 4 < 0 , roughly trac¬ 
ing the bulge component because they are along the minor axis of 
NGC 4697, all have generally older age and lower metallicity. The 
same holds also for points with larger | / 131 , although in this case the 
relations are less robust. 

Remarkably, the spatially resolved stellar population results 
perfectly agree (on a qualitative level at least) with the spatially re¬ 
solved kinematics results. Unfortunately none of the hydrodynam- 
ical simulations derive 2D detailed stellar population properties. In 
the future, it might be extremely useful to have predictions from 
simulations on the ages and metal abundances of the remnant of 
different mergers scenarios. 

In the S AURON ETGs of low to intermediate mass (cr = 100- 
160 kms“^), the region of increased metallicity is connected to 
a mild decrease of the [a/Fe], whereas in more massive galaxies 
(cr > 160 kms“^) the disk formation and increased metallicity is 
typically constrained to a central location involving only a minor 
fraction of the total mass of the galaxy and only in about 25 per 
cent of the cases also a central depression in [a/Fe] can be seen. 

In the case of NGC 4697, we find that the central region is 
slightly a—enhanced with respect to the bulge, whereas the disk 
shows roughly the same [a/Fe] as the bulge. A larger metallicity as¬ 
sociated with the fast rotating component, seems to point toward a 
formation scenario where fast rotator galaxies have undergone one 
or multiple periods of secondary star-formation at a later stage and 
over longer timescales and therefore resulted in an increased metal¬ 
licity and a decrease of the abundance ratio in a disk-like structure 
( [Kuntschner et al.[2010[ ). 

A younger and more metal-rich disk is also well consistent 
with the scenario suggested by the gas-rich merger (e.g. H-i-09), as 
the younger disk would arise when the gas, retaining its angular 
momentum during the merger, is fueled in the center of the rem¬ 
nant and eventually turns into stars. In Section [33] we have demon- 
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strated that this same scenario is also perfectly able to explain the 
kinematical substructures. 

We note here that the disk merging scenario does not neces¬ 
sarily imply major merging, as also merging with lower mass frac¬ 
tion can produce similar kinematical features. In fact, simulations 
in N-i-14 have shown that correlations similar to the ones we ob¬ 
serve in the hs and parameters of NGC 4697 can originate from 
two galaxy formation scenarios: a late gas-rich major merger or in 
situ star formation plus minor merging. 

In conclusion, the combination of the kinematical features and 
the stellar population properties seem to point to the scenario where 
some fast rotating and disky ellipticals can have formed from dis¬ 
sipative (wet) mergers of disk galaxies. 


5 DISCUSSION AND CONCLUSIONS 

In this paper, we have presented high S/N integral field kinemat¬ 
ics and stellar population combined analysis of the elliptical galaxy 
NGC 4697, an almost edge-on, intermediate luminosity (L*) iso¬ 
lated elliptical galaxy. We used a mosaic approach, the first of its 
kind with VIMOS IFU data, obtaining 2D maps of stellar kinemat¬ 
ics and absorption lines up to 0.7i?eff- The final effective field-of- 
view of ~ 80” X 60” allowed up to obtain one among the largest 
2D kinematical maps ever measured for a single early-type galaxy 
(in terms of arcsec) providing a unique view on the internal struc¬ 
ture of the system, as well as crucial morphological and dynamical 
information. 

Our main scientific findings are the following: 

i) We have measured the 2D velocity moments of the galaxy’s 
stellar population finding clear evidence for a rotation-supported 
structure along the major axis from the rotation velocity and the 
third Gauss-Hermite coefficient (representing the skewness). We 
have measured the “spin-parameter” (Ar) confirming the classifi¬ 
cation of this system as ‘fast-rotator’. 

ii) We have analyzed the correlations between the and 
parameters with V/a (Fig. [^, together with the Ar spatial pro¬ 
file, demonstrating how these characterize the different kinematical 
components (bulge and disk) of the galaxy. 

iii) We have directly compared these correlations with those pre¬ 
dicted by hydrodynamical simulations of [Hoffman et al.| ( |2009| ) 
from gas-rich disk-disk major mergers with a 30% gas fraction and 
of jNaab et al.| ( |2014) for gas-rich minor and major mergers (class A 
and B). The correlations we presented closely resemble the theoret¬ 
ical predictions and indicate that the mass assembly of NGC 4697 
is dominated by gas-rich minor mergers, gradual dissipation plus 
possibly a late gas-rich major merger. 

iv) A direct comparison with ID long-slit spectroscopy < [de| 
[Lorenzi et al.|2008) shows an overall good agreement between the 
different kinematics results along the major axis. 

v) We extracted ID spectra of the disk and the bulge of 
NGC 4697, using a two-components fit to a deep V-band image. 

vi) We studied the stellar populations of these two components 
using SSP models. We find indications that the stars in the disk are 
~ 2 — 3 Gyr younger and more metal rich that the stars in the bulge. 

vii) We linked the stellar population analysis and the internal kine¬ 
matics of the system by correlating age and metallicity of the stellar 
populations of the two components with the hs and the kine¬ 
matic parameters, finding further evidence for a significant in-situ 
formation of stars since z ~ 2 with additional gas-rich minor merg¬ 
ers (although we do not exclude the possibility of a late major- 
merger that left a previously rapidly rotating system unchanged). 


viii) From the two bluer IMF-sensitive optical indicators of S-i-14, 
which fall in the VIMOS wavelength range, we do not detect any 
difference in the IMF slope between the bulge and the disk. Both 
components are consistent with an IMF with a single power slope 
of x=2.0, consistent with a Kroupa IMF and slightly lighter than 
a Salpeter-like IMF. The inferred IMF slope is perfectly consis¬ 
tent with the one predicted from the IMF-a* rela tion for an ETG 


with cr* = 160 kms such as NGC 4697 ( 

jConroy & van 

Dokkum|20121|La Barbera et al. 

20131|Treu et al.||20101|Tortora, 

Romanowsky & Napolitano|2013 

ITortora et al.|2014b| S+14). 


The absolute novelty of this work resides: i) on the depth and 
width of the data acquired as we have very high S/N spectra over 
a an effective area of the sky of 80" x 60" (e.g. larger than MUSE 
FOV); ii) on the presentation of a detailed spatially resolved kine¬ 
matical and population 2D analysis within the central effective ra¬ 
dius, which has allowed us to cross-correlate LOSVD information 
(i.e. hs and h^) with relevant line indices (e.g. vs. H/3 and 
Mgb in Figure [^, showing that points with > 0, associate 
mainly to the later-formed disk, are younger and more metal rich, 
while points with /i 4 < 0, roughly tracing the bulge component be¬ 
cause aligned with the minor axis of NGC 4697, all have generally 
a smaller and smaller [M/H]. 

The emerging picture is that a later formed disk has been built 
from pre-enriched gas that formed stars after having dissipatively 
sunk into the center. This is well consistent with the scenario sug¬ 
gested by gas-rich disk merging simulations (e.g. H-i-09), as the 
younger disk would arise when the gas, retaining its angular mo¬ 
mentum during the merger, is fueled in the center of the remnant 
and eventually turns into stars. It is remarkable that within this sce¬ 
nario, photometry, stellar population features and kinematical fea¬ 
tures like the hs- and h^-V/a correlations are all qualitatively cor¬ 
rectly predicted. We stress again that H-i-09 only considers 1:1 disk- 
disk merger but that also merging with lower mass fraction (minor 
mergers) can produce similar kinematical features (e.g. N-i-14). Fur¬ 
thermore other mechanisms involving a rapid dissipative collapse 
of enriched gas (clumpy disks, violent relaxation, cold accretion), 
might have produced a similar galaxy component. 

This paper is intended as a pilot program for a large project 
that aims at the study of internal structure, mass distribution and 
stellar population of a sample of nearby L* ETGs. With medium- 
high signal-to-noise spatially resolved spectra, we will perform a 
detailed 2D study of the stellar population via line-index measure¬ 
ments and try to further test possible variation in the IME normal¬ 
ization. We will attempt to characterize the low-mass end of Initial 
Mass Eunction normalization in the bulge and in the disk compo¬ 
nents by measuring EWs of the IME-sensitive features defined in 
S-i-14. In order to achieve this result we plan to change the con¬ 
figuration of VIMOS in the new observations, choosing a redder 
filter (HR-red) to cover the proper wavelength rangj^ In addition 
to a direct measurement of the stellar mass-to-light ratios (M*/L), 
thanks to the 2D IE coverage, we will be able to study age and 
metallicity radial gradients for our systems. Moreover, where spa¬ 
tial gradients in the stellar population parameters are present, the 
assumption of a spatially-constant (M'^/L) clearly does not work. 


We are aware that moving to the red filter will cause the loss of some 
of the blue indices covered in this paper, making age determination more 
complex. We therefore plan to select ETGs where measurements of age and 
metallicity are available from the literature and mainly focus on measuring 
gravity-sensitive features to study the possible spatially variation of the IMF 
slope, with particular regard on bulge versus disk dominated zones. 
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Figure 10. Left panels: H/3 EWs versus the hs (top) and the /14 (bottom) parameter. Middle panels: Mgb EWs versus the hs (top) and the /14 (bottom) 
parameter. In all panels, points represents values obtained from the macro bins of the smoothed maps and show a correlation between LOSVD properties 
and stellar population features. Larger H/3 values roughly indicate younger ages, while larger Mgb values indicate higher total metallicity (although the 
well-known age-metallicity degeneracy makes things a bit more complicated than that; e.g. | Worthey | j 1 994) . Right panels: H/3 EWs versus the hs (top) and 
the /14 (bottom) parameter, color-coded by Mgb EW values. Points with larger, positive / 14 , aligned with the major axis (see Eig.|^ have on average higher 
and Mgb. This is true also for points with larger hs values (in absolute value) although the plot looks a bit noisier. This Eigure highlights the agreement 
between kinematics and stellar population analysis both confirming that stars aligned with the major axis (where the disk dominates) are younger and more 
metal-rich. 


Consequently, the 2D coverage, allowing a more accurate measure 
of the (M*/L), also permits a more accurate and complete dynami¬ 
cal analysis (e.g., |Tortora et al.|2010b|[20TT] ). 

A quantitative dynamical model to fit the whole 2D maps for 
the rotation and the velocity dispersion, is out of the purpose of this 
paper but will be performed in future work. Only with a complete 
2D kinematical map of a galaxy (up to at least as well as 

a robust dynamical model it will be possible to measure the DM 
fraction, trace its radial distribution as a function of the mass of 
the halo, and test plausible environment dependencies. In future 
work we therefore plan to perform a detailed study of the stellar 
content of ETGs to combine spatially resolved stellar information 
on the central region with planetary nebula data on the outskirts (up 
to ~ 5 ') with the final purpose of obtaining the complete rotation 
curve and distribution of dark matter (e.g. |Napolitano et aL|2009| 

I20TT] ). 

Moreover, thanks to a complete and detailed dynamical anal¬ 
ysis, we will have a completely independent measurement of the 


IMF slope that we will compare to the one obtained by stellar pop¬ 
ulation modeling. If both methods agree on the resulting IMF, this 
will be a strong test of the absence of hidden systematics. 


Finally another possible interesting follow-up will be study¬ 
ing the dynamics and the IMF normalization of this galaxy in the 
phenomenological framework of the Modified Newtonian Dynam- 


ics theory (MOND, Milgrom 

1983b|a| following the approach pre- 

sented in|Tortora et al.|2014b 
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